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Abstract 

We introduce a new matter action principle, with a wide range of applicability, for the Vlasov 
equation in terms of a conjugate pair of functions. Here we apply this action principle to the 
study of matter in Bianchi cosmological models in general relativity. The Bianchi models are 
spatially- homogeneous solutions to the Einstein field equations, classified by the three-dimensional 
Lie algebra that describes the symmetry group of the model. The Einstein equations for these 
models reduce to a set of coupled ordinary differential equations. The class A Bianchi models 
admit a Hamiltonian formulation in which the components of the metric tensor and their time 
derivatives yield the canonical coordinates. The evolution of anisotropy in the vacuum Bianchi 
models is determined by a potential due to the curvature of the model, according to its symmetry. 
For illustrative purposes, we examine the evolution of anisotropy in models with Vlasov matter. 
The Vlasov content is further simplified by the assumption of cold, counter-streaming matter, a 
kind of matter that is far from thermal equilibrium and is not describable by an ordinary fluid 
model nor other more simplistic matter models. Qualitative differences and similarities are found 
in the dynamics of certain vacuum class A Bianchi models and Bianchi Type I models with cold, 
counter-streaming Vlasov matter potentials analogous to the curvature potentials of corresponding 
vacuum models. 

PACS numbers: 04.20.Fy, 98.80.Jk, 52.25.Dg, ll.10.Ef 



2 



I. INTRODUCTION 



According to observations, our Universe is both spatially homogeneous and isotropic on 
large scales [TJ [2]. However, in this paper a more general case is considered in which the 
universe is taken to be spatially homogeneous but anisotropic [3j 0]. We further assume 
that our models contain matter governed by the Vlasov equation [HH7], which opens the 
possibility that it can be anisotropic in momentum space as well. Spatially-homogeneous 
cosmological models, which are also known as Bianchi cosmologies [3] , are classified into nine 
types based on a standard classification of three-dimensional Lie algebras [S] that determine 
the symmetry of the model. 

The presence of matter influences the dynamics of the model. In this paper we study 
Vlasov matter: collisionless particles interacting only through their mutual gravitational 
effect. We introduce a new action principle for such matter, applicable in a wide variety of 
contexts. Hamiltonian methods have long had important roles in general relativity [91 [10] . 
In a model that is filled with Vlasov matter in a spatially-homogeneous way, the presence 
of matter appears as an additional potential term in the Hamiltonian [TT] . These additional 
potentials are not as steep as the curvature potentials of the vacuum Bianchi cosmologies, so 
the evolution of anisotropy is always under the influence of the matter potential, in contrast 
to the vacuum cases [TT]. In this paper, vacuum Bianchi cosmologies and Type I models 
with Vlasov matter are compared in order to clarify the distinction between similarly shaped 
curvature and matter potentials. 

There are two main reasons for studying Bianchi cosmologies. First, universes that 
are spatially homogeneous but anisotropic are the simplest generalization of the spatially- 
homogeneous and isotropic universes, because there are no models that are everywhere 
isotropic and spatially inhomogeneous [3] . From the theoretical point of view, it is of inter- 
est to see how the dynamics of a model universe is changed by this first step away from the 
spatially isotropic cosmologies. Second, due to the spatial homogeneity, the Einstein equa- 
tions are reduced to a set of coupled ordinary differential equations, which can be viewed as 
a finite-dimensional dynamical system. 

A reason for studying matter described by the Vlasov equation is that this model allows 
for phase space degrees of freedom. Consequently, with it the dynamics of nonthermalized 
matter can be explored. The reduction to a set of ordinary differential equations is also 
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possible with this kind of matter description within the Bianchi cosmology context. We 
believe that physically meaningful cosmological models should have realistic matter models. 
Although both the Bianchi cosmologies and the Vlasov matter we consider are simplistic, 
they are less so than the common practice of merely describing matter dynamics by a perfect 
fluid. We thus agree with several others in this goal (see [12]). 

A spatially-homogeneous spacetime is foliated by spacelike hypersurfaces, on which the 
metric tensor g is invariant under the transitive action of an isometry group. If the group 
is a simply transitive Lie group (and thus three-dimensional), symmetry is characterized by 
three linearly- independent, spatial Killing vector fields {$,(i)}^ =1 . Their constant structure 
coefficients C \ • determine the nature of the symmetry: 

where the Einstein summation convention is used, as it will be throughout the paper. The 
three-dimensional Lie groups (actually their Lie algebras) are classified into nine Bianchi 
types. This classification of the symmetry of the spacelike hypersurfaces defines the Bianchi 
cosmological models, which were further subclassified as class A and class B models by 
Ellis and MacCallum [13]. The models of interest in this paper, the class A models, are 
characterized by C k ik = 0. The structure coefficients used in this paper for the class A 
models are in common standard forms that are summarized in Table IB 

TABLE I: The structure coefficients (structure constants) of class A Bianchi models [4j. These 
coefficients are given in standard form and of course may vary under linear transformations of the 
basis vectors. 



Type Structure Constants 

I C% = 

II C\ 3 = -C\ 2 = 1 (other: 0) 

VIo C\ 3 = -C\ 2 = 1, C7 2 13 = -C7 2 31 = 1 (other: 0) 

VIIo C\ 3 = -C\ 2 = -1, C7 2 13 = -C7 2 31 = 1 (other: 0) 

VIII C\ 3 = -C\ 2 = -1, C7 2 31 = -C7 2 13 = 1, C\ 2 = -C 3 21 = 1 (other: 0) 

IX C% = e ijk (the Levi-Civita symbol) 



When a class A Bianchi-type symmetry is imposed on the vacuum Einstein equations, 
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the equations become a coupled set of ordinary differential equations that can be written 
in Hamiltonian form. The potential terms arise from the curvature of the model which is 
dictated by the symmetry. The curvature potential terms in the various Hamiltonians (see 
Table [IT]) result in qualitatively distinct evolutions of anisotropy in the vacuum models [I]. In 
order to discuss the influences of the curvature potentials, it is necessary to understand the 
dynamics of the model universe in the absence of the curvature potential terms, namely the 
vacuum Bianchi Type I model. For this model the corresponding Hamiltonian is integrable, 
and the explicit solutions are known as the Kasner solutions. In the other Bianchi models the 
presence of a curvature potential affects the dynamics. The Kasner solutions are important 
because the dynamics of the other models can be approximated as a series of Kasner solutions 
with different Kasner parameters pE]. The curvature potentials are so steep that they can be 
approximated as moving potential walls within which the point that represents the state of 
the universe moves as an approximate Kasner solution whose parameters change when the 
universe point bounces off the walls. This approximation is known as qualitative cosmology 

HEMS]. 

The paper is organized as follows: In the next section, the Hamiltonian approach to 
Vlasov matter, based on a pair of conjugate potential functions, is presented in a gen- 
eral context. This variational principle, although based on earlier works [TT14T9] . is new. 
Spatially-homogeneous universes are discussed in a general way, and the Vlasov Hamilto- 
nian is given in a form compatible with the symmetry. In the ensuing section, the evolution 
of anisotropy in the presence of Vlasov matter is analyzed by comparing and contrasting cer- 
tain vacuum Bianchi models with those having Vlasov matter in a Type I model. (A cold, 
counter-streaming distribution function that supports the spatial symmetry is assumed.) 
The last section is devoted to conclusions and suggestions for further research. 

II. FORMULATION 

A. Variational Principle and Derivation of Hamiltonian 

In Vlasov theory [T9H22] . matter is modeled by a phase-space distribution function de- 
noted by F{x^,p v ), where the {x^} are the positions of the particles and the {p v } are their 
4-momenta. The systems that are studied here consist of particles of the same mass, so that 
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a mass-shell constraint will be imposed: 

g^PvPiy = ~ m2 and p° > 0. (2) 

For the purposes of this exposition, we set a gauge condition on the spacetime metric g^ v 
such that goo = — 1 and g i — for i — 1, 2, 3 (we later generalize this gauge). The invariant 
volume element in 4-momentum space is reduced to the following volume element on the 
3-momentum mass-shell, with components pf. 
d 3 p 



with p° := a/ m 2 + g^PiPj and g := det {g^l- (3) 



Because of the gauge choice and the mass-shell constraint, it is convenient to define an 
on-shell distribution function, 

f{x^p i ):=F{x^p i M^,Pi)), (4) 
where po = —p° is the momentum defined in equation (J3J). The notation p^ will be used to 



indicate the three variables pi = pi along with the functional form of po = —p°. 

Since matter is assumed to be collisionless, F is constant along geodesies, and so F and 
/ are governed by the relevant Vlasov equation. The off-shell Vlasov equation is simply the 
geodesic equation for the particle paths (the momentum one-form p is the mass m times 
the particle velocity): 

The on-shell Vlasov equation is obtained by using equation (J2|. It can be expressed as 
where {•, -}3 is the three-dimensional Poisson bracket: 

tf -01 1 f-i 9 f 1 a/3 - - 9f\ 

The metric tensor g in the Vlasov equation is determined by the Einstein equations 
(where the coupling constant is taken to be unity): 

Rfiu - ^Rg^ + A fiV = T ^ '■= J J^-Q p^Pvf^iPi) ■ (8) 

The action functional that gives equation Q is 

S = Sffilbcrt + + Setter = J d*X \/~gR - 2 A J d 4 X y/^g + Matter- (9) 



The Hilbert action Snubert gives the Einstein tensor and Sa gives the cosmological term 
upon variation with respect to g^ v and when necessary integrating by parts and ignoring 
the boundary terms [TJ |2]: 

1 ^Sffilbcrt 1 1 SSa 

j=§-6j*T- = ~ 2 R9 ^ and 7^gJg^ = Kg ^ (10) 



Matter must satisfy 

1 <5*Smatter , , , 



V^9 Sg^ 

An action using Matter must also yield the Vlasov equation. 

In order to put the Vlasov equation into a variational form, we express the distribution 
function by a conjugate pair of phase-space functions (A4,Af), using the four-dimensional 
Poisson bracket {•, -} 4 [TTHT9] : 

F(^^) = { MAr }4 :=—--— — . (12) 

M.{x^,p v ) and N{x^ip v ) are each required to satisfy Vlasov equations, which are compactly 
written as 

The Vlasov equation for F{x^,p u ) [23J, 

{f, \g^p^} A ={{M^U ^PmP^} 4 = o , (14) 

is derived by making use of the Jacobi identity [T8l [19]. The on-shell restriction can be 
achieved by defining 

fx(x tM ,p i ) := Mix^TPupoix^iPi)), v{x»,pi) :^Af(x ti ,p i ;po(x^,p i )). (15) 

Therefore, the on-shell distribution function is expressed by the on-shell pair: 

f(x^ Pi ) = {^u} 3 . (16) 

The equations that fi and v have to satisfy are also derived by imposing the mass-shell 



condition on equations (13), 



d I fJ> \ \ { P 
dt , „ , I \ v 



f) =0, (17) 



from which the on-shell Vlasov equation for f(x^,pi) is derived: 

-{^z/} 3 +{{/i,z/}3,p } 3 = 0. (18) 

Again, the derivation rests on the Jacobi identity. 

There is a gauge group associated with the introduction of the pairs (Ai,J\f) or (/x, v). For 
example, a transformation of the form (p,(fi, v\ u(fi, v)) that satisfies the Jacobian condition 
<9(/i, z/)/<9(/i, z/) = 1, ensures that {/2,z/} 3 = {fi, A full discussion of the gauge group 
associated with this 'potential' decomposition, including both dependent and independent 
variables, will be presented in a subsequent paper. 

Now that the distribution function is expressed by the introduction of a conjugate pair of 
phase space functions, the matter action for the Vlasov-Einstein system can be written down. 
The action principle in terms of these variables is canonical, that is the evolution of those 
variables is manifestly determined by a canonical Poisson bracket. We check the variational 
formulations in both off-shell and on-shell forms and describe how they are naturally related 
to each other. 

If there is no mass- shell constraint, the matter action is simply 

S oS [M,Af,g^} = - f d A x d 4 pg p(7 p pP(T {M,Af} 4 , (19) 



and the functional derivatives are 

1 ^off f d 4 p 



-9 8g 



Hi/ 



p„p„{M,AfU = -T^, (20) 



and 



0. 



6S&/6M = o] | (p^-lg^p.p^)^ , 

5S oS /5M = 0j = 0. 



Since we consider a system of equal mass particles, (19) has to be put on-shell. This can 
be done by restricting M. and M to be on-shell and integrating (19) over p Q , either before 
or after the variation. (In other words the on-shell constraint and the variation commute.) 
The integration over po is best achieved with the use of the function [2] 



S + (-g^p, Pu - m 2 ) := 2e(p°)6(-g^ Pll p v - m 2 ) = ^5(p - p ). (22) 

However, care should be taken with the independent dynamical variables. If variations 
are done before the on-shell constraint, then the independent dynamical variables are of 
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course g^ u ,M, and A/", but if variations are done after the on-shell constraint, they are 
g^,yU, and v. If the variational derivatives are taken before the on-shell constraint, the 
functional derivatives are found to be 



5S 



and 



5S/5M 
5S/5N 



d 4 p 
\f~9 



P»Pv {M,Af} 4 5 + (on shell) = —T f 



(23) 



-.F(x*,p v ) 



N\ 5+ (on shell) 
m) 5+ (on shell) 



0. 
0. 



(24) 



The on-shell delta functions in (24) can be eliminated by integrating them with respect to 



Po, which results in (17). The calculation of (23) is not as straightforward as (20), because 



the g^ u in 5+ (on shell) is subject to the variation. The calculation is lengthy, but all the 



extra terms vanish to give ( 23 ) [TT] . 

On the other hand, the mass-shell constraint can be imposed before variations. For this 
case, the action is given by 

dfj, 



S on [pi, u, <T] = 2 d*x d s p !/-£ - f{fl, u} 



dt 



whose functional derivatives are 

1 5S on 



-9 8g 



d 3 p 



-gp 



(25) 



(26) 



and 



SS on /5u 



dv/dt + {v,fy. 



(27) 



! = 0, 
dfi/dt + {^p°} 3 = 0. 
Note that the integrand of (25) is reminiscent of the phase space action, u pq — H" , since 
/i and v are canonically conjugate to each other and p° 



— p is the energy of a particle 



with u} 3 being the particle distribution on phase space. Consequently (27) have the 
Hamiltonian form: 

dfi _ 5H ^ dv _SH 
dt 5v dt 5/U 



(28) 



B. Vlasov Matter in a Spatially-Homogeneous Spacetime 



Since the spacetime structure and the matter configuration are related by the Einstein 
equations, the matter term is subject to symmetry conditions. If some spacetime symmetry 
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is assumed, the form of a distribution function is restricted so that its corresponding stress- 
energy tensor has the same symmetry as the spacetime. It is sufficient, in the presence of a 
Killing vector field that [25] 

OF dF 

^-^>= a (29) 

The quantity 

Y := p(£) (30) 



is a constant along the geodesic particle paths. If the Killing vector is spatial, equation (29) 
becomes 

{/,n3 = {/,p(*)}s = o. (31) 

Under spatial homogeneity, it is natural to use a left-invariant vector basis {X( M )} and its 
dual basis {a^'} [I]. The vector basis is chosen so that the spatial {X(j)} are tangent to 
the spatially- homogeneous sections, and {X(o)} is perpendicular to them. The vector basis 
satisfies the commutation relations 

[X ( o),X w ] = 0, (32) 
[X (i) ,X (i) ] = -C"yX w . (33) 

Similarly, the dual basis (with cr^ = dt) satisfies 

dcr(°) = (34) 
dcr« = ^ffWAffW. (35) 

The components of the metric tensor in this basis are functions only of t: 

9i*u = S ( X ( M ), X H ) = g^(t). (36) 

Momentum is expressed as 

h ll = p(X 0t) ) i.e.,p = h^\ (37) 
Note, however, that the hp are not canonical coordinates in general, because 

{h i ,h j h = C l ij h ll (38) 
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while canonical coordinates would satisfy {pi,Pj}3 = 0. With the momentum coordinates hi 
and the time variable t, the spatially-homogeneous distribution function has the form, 



f = f(t,K), (39) 
so that / does not depend explicitly on the spatial coordinates. The first term of (25) gives 



zero upon variation with respect to g^, and {//, z/} 3 can be replaced by /. The Vlasov 
equation can therefore be rewritten, reflecting the symmetry, as follows: 

df df „ d r dh° 



c)t + dt Cdabkd dh = °' ^ h ° = = V'" 2 • U" h > h J- ( 4 0) 

It is immediately seen that in Type I (where C c ab = 0), the solution of the Vlasov equation 
is independent of time: / = f(hx, h 2 , h 3 ). 

Based on the matter action we derived in the previous sub-section, the Hamiltonians 
for the class A Bianchi cosmologies with Vlasov matter are obtained. The calculation is 
simplified by restricting to the diagonal case: The positive-definite spatial metric g^ will be 
taken to be diagonal in a basis that is invariant under the symmetry transformation described 
by the Lie group of the particular Bianchi model [26]. This is not a gauge choice; it imposes 
physical limitations on the model (see p31 E3 EH]) not only because the standard form for 
the structure coefficients of the Lie Algebra are used, but also because of restrictions which 
must be placed on the particle momenta (as described below). At this point it is necessary 
to keep all diagonal components of the metric as variables (rather than setting g 00 = — 1 as 
we had done earlier) in order to formulate explicit Hamiltonians: 

9lxv = g(X M , X H ) = diag {-N 2 (t), A 2 (t), B 2 (t), C 2 (t)) . (41) 

As we said, this restriction to the diagonal case, when the standard forms from Table |T] are 
used for the C" jfc , is really a restriction in a sense that some physical effects are then not 
allowed. The T^ v based on the distribution function must also be diagonal; for example, 
it is sufficient that the distribution function be an even function of hi,h,2, and h^. After 
obtaining the Hamiltonian, we can specify the form of N, by choosing an appropriate time 
coordinate. Thus, the Hamiltonian systems we obtain in this section are three-degree-of- 
freedom systems. Note that the components of the metric tensor depend only on t due to 
the use of the invariant basis. 
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The action functional is 



Smibert + S A + S mattcr = J dt ~ (ABC + ABC + ABCj + (terms with C l jk ) 

/f f I L~2 h/2 i 2 

dtNABC-2jdt J d 3 hNf(t,h k )^m 2 + ^ + ^ + -^, (42) 

where ' is d/dt and where "terms with C appear as a curvature potential term in the 
Hamiltonian. By performing the Legendre transform on the Lagrangian, the Hamiltonian 
is obtained. Since there is no N, N is treated as a Lagrangian multiplier, and the variation 
with respect to N gives the Hamiltonian constraint: 

= H:= (A 2 tt 2 a + B 2 tt 2 b + C 2 n 2 c - 2ABtt a tt b - 2BCn B n c - 2CA<k c >k a ) 



// U 2 U 2 U 2 

d 3 hf(t,h k )^jm 2 + ^ + ^ + -^. (43) 

It is often convenient to use the variables a(t), /3±(t), which are defined by: 

In terms of these coordinates, the curvature potentials are as shown in Table [TT} Note that 
the dependence on a is always of the form V c = e a v c ((3±). The coordinates of (44) are 
called the Misner parameterization [29] . The physical meaning of these dynamical variables 
is clear: e a is the universe scale factor, and the /3± characterize spatial anisotropy. In terms 
of the Misner parameterization, the class A Hamiltonians have the following form: 

H = N (J^e~ 3a (-7T* + ir} + + tt|_) + 2Ae 3a + V C + V m j . (45) 

The equipotential curves for V c in the /3 + — (3- -plane are depicted in Figure [TJ 

The term that comes from the presence of Vlasov matter is also a potential, which is 



denoted by V m in equation (45). The matter potential is 

V m = 2e~ a J d 3 hf^m 2 e 2a + (h 2 e ~ 2 ^- 2 ^- + h 2 e - 2 ?++ 2 ^- + h 2 e 4 "+) . (46) 
In massless cases, the dependence on a can be separated from that on (3± since 

V m = e~ a v m (P±) (47) 

where 

V m = 2 J d 3 h fx/h 2 e -20+-2>/3/3- + h 2 e -2p + +2V3/3- + fc2 e 4/3+ _ ( 4 g) 
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TABLE II: The curvature potentials. These potentials are calculated using a diagonal metric in 
a group-invariant basis: = diag(— N 2 (t), A 2 (t), B 2 (t), C 2 (t)) . The a(t),/3±(i) are defined by 
A = e a+A++v30- } b — e a+/3+-V3f3- ^ anc j q _ g Q-2^ + _ (Note: f or calculations using a general form 

of the metric, see |30j.) 



Type 



V c = e a v c ((3±) 



I 

II 
VIo 
VIIo 
VIII 

IX 



V" - 2BC 

tfa = aZSc + B2 f = 2e Q e 4 ' 3 + cosh 2 (2 v / 3 / S_ 
VV// = i A2 -B 2 f = 2e a e^+ sinh 2 (2V3/3_ 



t/ _ 1 / A 3 , B 3 i C 3 \ /AB BC CA-i 
™ ~ 2 [BO + CA + AB ~ [-C A B~) 



= e a (2e- 2l3 + cosh(2V3/3_) - e 4/3 + + e 4/3 + cosh(4 v / 3/3_) + ^e" 8/3 +) 

t/ _ 1 / A 3 i B 3 , C 3 \ MB , BC , CAN 

K ^ - 2 (,bc + ca + as) ~ Kir + "a" + -b) 

= e a (-2e" 2/3 + cosh(2 v / 3/3-) - e 4/3 + + e 4/3 + cosh(4^/3/3_) + 5e" 8/3 +) 



Type Vll 




FIG. 1: Equipotential curves for class A curvature potentials. 

III. RESULTS: VLASOV EFFECTS 

In this section, we give a description of how a cosmological model containing Vlasov 
matter may mimic or may differ from an apparently similar vacuum model. The results will 
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be described, but the details of the calculations — both analytical and numerical — will be 
left for a subsequent paper. The general behavior in the sense of qualitative cosmology will 
be described, where the effects of the potentials are either negligible or act as hard walls. 
Keep in mind that these results were obtained for the case where A = and m = 0. 



From the Hamiltonian (45), we see that 
dH 



N 



a 



Pi 



Tin 



dn a 
dH 

dH 

da 



12 

N 



e 7T r 



12 



e vr /3± , 



N 



-3a 



-7T, 



7T 



2 

0+ 



IT, 



) 



d_ 
da 



(V m + V c ) 



(49) 
(50) 

(51) 



Presuming that a is monotonic, a is used as the "time" variable. In that case the dynamics 
of the model will be governed by the behavior of (3±{a). The motion of the universe point 
(a point in the (3 + — /3_-plane) is with speed w, using a as the "time variable" : 

d/3 + \ I d/3- 



w 



da 



da 



The Hamiltonian of equation (45) is used in the Hamiltonian constraint (43) to find 

2 -l-2Ae 3a (V c + V m )/7r 2 a . 



w 



(52) 



(53) 



When the potentials may be neglected, w = 1; this motion corresponds to the Kasner 
model. In this sense the models may be thought of as a series of Kasner epochs which tran- 
sition from one Kasner state to another when the universe point interacts with a potential 
wall. Note also that when the potentials may be neglected, 7T a is proportional to H and 
therefore vanishes. 

The vacuum Bianchi models of Types II, VIo, and IX will be compared with Type I 
models containing Vlasov matter. (See also [3 X J . ) Since the Type I models do not have 
a curvature potential, this procedure emphasizes the effects of the matter potential. From 
Figure [TJ we see that the qualitative effects of the curvature potentials in the vacuum models 
II, VIo, and IX are characterized by one wall, two walls, and three walls, respectively. The 
equipotentials of Figure [T] define these walls, which move in time, and this movement will 
be described in more detail later. 

For the matter potential, the distribution function / (which in the Type I case is inde- 
pendent of time) is chosen to be that of cold, counter-streaming matter: 
K 



5(h x -a) + 5(hi + a) 5(h 2 - b) + 5(h 2 + b) 5(h 3 - c) + 5(h 3 + c) 



(54) 
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An example of the potential functions fi, v of equation (16) may be easily calculated in this 
case, since in Type I, the hi are canonical variables: 



x 1 
K 



(55) 

H(h x - a) + H(h + a)] \s(h 2 -b) + 5(h 2 + 6)1 \s(h 3 - c) + S(h 3 + c)l , (56) 



where H is the Heaviside or step function. Of course, other forms for /x, v which yield f(hi) 
are possible, and the gauge group of these possibilities will be presented in a subsequent 
paper. 

V m for this model is 



V m = 2Ke- a Vm 2 e 2a + a 2 e-^+- 2 ^- + fe2 e -2/3 + +2 v / 3/3_ + c 2 e 4/?+ 



(57) 



This distribution function is for eight streams of particles, counter-streaming by pairs. The 



choice (54) for / is motivated by the cold, two-stream instability of plasma physics, where 



a large amount of free energy is stored in an equilibrium state that is released by linear 
instability. Studies of counter-streaming matter in the Newtonian gravitational case show 
that similar instabilities exist, and even explosive instabilities may be expected [321 33] • 
Although this particular choice is still simplistic, it is consistent with the Bianchi symmetry 
and, importantly, it provides a tractable example of the kind of effects to be expected when 
matter is nonthermal and anisotropic in momentum space. We note that although individual 
beams are cold, this model does have a nonzero stress-energy tensor, which, because this 
distribution function is even, is diagonal. 

If a = b = 0, c > 0, the model is called one- wall Vlasov. If a = b > 0, c = 0, it is called 
two-wall Vlasov. Finally, if a = b = c > 0, it is called three-wall Vlasov. These are the only 
models considered for inclusion in this paper, with m = just for the sake of illustration. 
The forms of the matter potentials, labeled Vi, V2, V3, are illustrated in Figure [2} 



One-Wall Vlasov 



Two-Wall Vlasov 



Three-Wall Vlasov 



-2 n 2 





FIG. 2: Equipotential curves for the Vlasov-matter potential with cold counter-streaming matter. 
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These potentials are given analytically by the following: 
One- wall Vlasov: a = b = 0, c > with m = 

Vi = 2Ke- a ce w+ . (58) 

Two- wall Vlasov: a = b ^ 0, c = with m = 



V 2 = 2V2Kae- a e~P + \J cosh(2 VZ/3- ) . (59) 
Three-wall Vlasov: a = b = c ^ with m = 



l/ 3 = 2Kae- a \j2e- 2 P+ cosh(2v / 3/3-) + e 4 ^+ . (60) 

These models of Vlasov matter in a Type I cosmology will be compared, respectively, to 
vacuum Bianchi models of Types II, VIo, and IX (compare Figures [T] and [2| . The curvature 
potentials for the vacuum models are labeled Vn, Vyi, and Vix] see Table [TTj The curvature 
potential for a Type I cosmology vanishes. Therefore in all cases when the universe point 
has an energy substantially above the value of the potential, the model acts like a Kasner 
model. The similarities and differences in the three situations described depend on whether 
the universe point is strongly affected by the potential, in which case the universe point is 
said to bounce off a potential wall or be influenced by it. 

The first pair is vacuum Type II compared to one- wall Vlasov (in Type I). The one- 



wall effective potential V e s, from equation (53), is proportional to e 2a+2l3+ , so that the 
equipotential walls move with speed 

w x = 1. (61) 
In contrast, the effective potential in the vacuum Type II case, from Table [TTJ is proportional 

tO e 4«+4(/3++v / 3/3-) ; go that the gpeed of the waU ig 

mi = \ (62) 

Consequently the universe point in the vacuum Type II case will (in many cases) strike a 
potential wall and bounce (once) in a very well-defined manner and then move away from 
the wall. Qualitative cosmology is useful in this case. In contrast, the universe point in 
the one-wall Vlasov model and a potential wall both have the same speeds. If the universe 
point is not moving exactly in the opposite direction of a wall, it will bounce. After the 
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bounce, the universe point cannot escape the influence of the potential, since equation (53) 
shows that its velocity will in general be less than 1. Therefore qualitative cosmology is not 
a good approximation. (Note, the one-wall case corresponds to an invariant boundary in 
the analysis of [3l].) 

The two-wall Vlasov model and the Type VIq model also have a qualitatively similar 



potential wall structure. The effective potentials for these two cases, from equation (53), is 
e 3a times the potentials listed in Table [TT] or equation (59): 



V$ = 2e 4a e 4/3 + cosh 2 (2>/3/?_) (63) 



V 2 eff = fce^e^+Wcosh (2^/3-) . (64) 
The speeds of the vertices of the equipotentials (that is, when /3- = 0) are: 

w v v f = 1 (65) 

w v 2 ert = 2. (66) 

The speeds of the walls, however, are governed by their slopes far from the vertices. At large 
/?_ the effective potentials are approximated by 

V?/ « e 4a e 4/3++4 ^- (67) 

^//~^ e 2a e -^ + +^_ (6g) 

The speeds of the walls, therefore are 

1 . , 

Wvi = r (69) 

w 2 = 1 (70) 

(as in the previous case). Therefore the same comments about the universe point hitting 
the potential walls apply: Qualitative cosmology is a good approximation for the Type VIo 
case but not for the two-wall Vlasov case. (Note, as with the one wall case, the two-wall 
case corresponds to an invariant boundary in the analysis of |31j.) 

The three-wall Vlasov and the Type IX vacuum Bianchi models are somewhat more 
interesting. The Type IX model eventually collapses. Note also that the potential walls have 
channels (see Figure [T]), which together with the steepness of the walls produce what Misner 
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calls the Mixmaster Model [29]. The three- wall Vlasov model has an effective potential of 
e 3a times the potential of equation (60): 



V 3 eff = ke 2a ^2e-^+ cosh (2^/3-) + e^+. (71) 

The speed of each vertex (for example when /3_ = 0, along the negative /3+-axis) is 

w v 3 ert = 2. (72) 

The speed of each wall is found by looking at the vertical wall when /3_ = 0, along the 
positive /3 + -axis: 

^3 = 1 (73) 

Again it is seen that the details of the motion of the universe point cannot be approximated 
by qualitative cosmology. 

The details of the motions in the three Vlasov models as contrasted with the correspond- 
ing vacuum models, in spite of the apparent similarity of the potentials, will be left to a 
subsequent paper. However, in the two-wall case the speed of the vertices of the equipoten- 
tials shows that partial isotropization takes place (namely — > 0). The three- wall Vlasov 
model has complete isotropization: The universe point goes to the origin. 

To examine this isotropization more closely, we first look at the stress-energy tensor, 
equation ([8]), for the three- wall case in a Type I model. In this case the pi are the same as 
the hi, and in our gauge, sf—g = NABC. Moreover, the on-shell condition implies 



P° 

y N 



IfVri* + (Pi/A) 2 + (P2/BY + (p 3 /Cy. (74) 



The stress-energy tensor is diagonal, with components: 



T 00 = iV 2 J^IT (75) 

T ^ A% ^cW^ (76) 

T22 = b2 abciwi (77) 

K c 2 

T33 = c2 abciic^ (78) 



where the constants K, a, b, c are defined in equation (54) and where 

IT := Vm 2 + {a/A) 2 + (b/B) 2 + {c/C) 2 . (79) 
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In the limiting case of A = B = C (when /3± = 0), and if it is approximately true that 
a = b = c, then the stress-energy tensor becomes that of a perfect fluid with W the energy 
density and V the pressure, given by 



Their ratio is the equation of state: 

V 1 

W = 3 + (mA/a) 2 ' (81) 
If m 7^ 0, the pressure goes to zero at large A (so-called "dust"), and if m = 0, the factor 
1/3 is that of a fluid of massless particles. This is consistent with results of j3TJ G3] that 
indicate that models of a class, which includes our three-wall model, will isotropize in the 
asymptotic future. 

As for the action integral, the h integration can be done to get, in our diagonal Type I 
model, 

^matter = -2K J NYidt = -2K J Nyjw? + (a/A) 2 + (b/B) 2 + (c/C) 2 dt 

= -2K J -^=y/m 2 + a 2 g u + b 2 g 22 + c 2 g 33 dt. (82) 

The variations with respect to g^ v give the stress-energy tensor components displayed before. 



IV. CONCLUSION 



In this work, the matter action for the Vlasov-Einstein system for equal mass particles 
is constructed. As an illustration of this action, we found a new kind of potential V m for 
Vlasov matter in the anisotropy plane of spatially-homogeneous models (those of Bianchi 
class A with diagonal metric). 

Then, in order to investigate the similarities and the differences between the Vlasov- 
matter potential and the curvature potential, the Type I universe with cold, counter- 
streaming matter and the corresponding vacuum Bianchi models were compared on the 
basis of the shape of the potentials. Both kinds of potentials, the curvature potentials V c 
which arise from the requirement of spatial homogeneity in a vacuum model and the V m , were 
classified by the number of approximately straight equipotential sides. The V c in models of 
Type II has one; in Types VI and VII the V c have two; and in Type VIII and IX the V c 
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have three. Their counterparts, one-wall, two-wall, and three-wall Vlasov-matter potentials 
were constructed by selecting a cold, counter-streaming matter distribution function in a 
Type I model. The V c of a vacuum Type II model was compared with the one-wall Vlasov 
model. The V c of a vacuum Type VIo model was compared with the two-wall Vlasov model. 
The V c of a vacuum Type IX model was compared with the three-wall Vlasov model, since 
they share triangular symmetry, in spite of the absence of channels in the Vlasov case. 

Aside from questions of recollapse and direction of increase or decrease in anisotropy, 
it was seen that the V c and V m had different effects because of the relative speeds of the 
walls and the universe point described by /3±(a). The dynamics of the vacuum Bianchi 
cosmologies can be approximated as a series of Kasner eras where the point representing the 
metric in the (3± space moves freely except for bounces off the potential walls. Because V m 
is not as sharp as V c , and because the speed of its walls was comparable to the speed of the 
universe point, it would be difficult to separate out a potential- free region: Consequently, 
contact with V m could not be characterized as a transition from a Kasner state to another 
Kasner state. In a subsequent paper, details of the effects of Vlasov matter on the evolution 
of anisotropy will be described both for our present case of cold, counter-streaming matter 
and matter which may be considered warm. 

Early work on this subject was done by Misner [35] and Matzner [36]. Their idea was to 
ask whether a kinetic theory of matter would cause anisotropy to decay faster than it would in 
a model with perfect fluid. They took a general approach, in which the distribution function 
had a thermal structure because their massless particles (which they called neutrinos) were in 
thermal equilibrium with other matter until a particular time when the temperature dropped 
below a critical value. Then there was a period of transition followed by the present epoch, 
when the neutrinos are essentially collisionless. Misner looked at Type I models. Matzner 
looked at Type IX and Type V models, with and without rotation. They were able to 
get a general idea of a potential due to the particles. This potential was a triangular one, 
possibly with fewer than three walls, depending on the form of the distribution function. 
Our methodology differs from theirs. Our distribution function in our main example also 
differs, but our general results are in accord with theirs. 

Rendall |34j has also studied the effects of Vlasov matter in a Bianchi Type I model. 
His results cover a more general set of distribution functions than we considered in our 
particular example. In so far as our results can be compared, they agree. Rendall did 
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not use a Hamiltonian formalism, and we feel that subsequent studies will benefit from our 
approach. 

In conclusion, we list several possible areas for future work. Examples of minor projects 
are as follows: First of all, the distribution functions of Ehlers type were used in this work, 
and the dynamics of the more general distribution functions given in equation (3.25) of 
Chapter 3 of [11] should be treated. It would be of interest to perform numerical compu- 
tations for these more general f(t,hk)- Second, our system was simplified by assuming a 



diagonal metric (41), and it could be generalized to include off-diagonal elements and possi- 
bly even rotating models. Third, our examples focused on massless cases for simplicity, and 
the massive case should be investigated. 

There are at least two major projects. First, we propose to study further the problem of 
a Bianchi Type IX model with non-zero cosmological constant and Vlasov matter. This case 
has Einstein's static universe as an unstable fixed point. Whether this point is a chaotic 
scatterer or not for dust filled models was discussed in [3"TH4"T] . which has no additional 
anisotropic potential. The inclusion of V m for Vlasov matter might change the behavior 
of the scatterer, and introduce possible chaotic effects. Second, we propose the problem 
of the Type I universe with the Vlasov-matter potential V m in the presence of a negative 
cosmological constant. In this case, the universe eventually re-collapses and the potential 
walls recede, just as for the vacuum Type IX case. However, unlike what happens in the 
vacuum Type IX case, the universe point cannot catch up with the retreating potential 
walls of V m , since they can move with a speed comparable to that of the universe point. 
Therefore, with a negative cosmological constant, the system can be characterized as a 
transition from an initial Kasner state to a final Kasner state with the intermediate dynamics 
being complicated: that is, the system can be viewed as a scattering problem. 
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